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Abstract 

In order to study further the long-range correlations ("ridge") observed recently in p+Pb collisions 
at ^/snn = 5.02 TeV, the second-order azimuthal anisotropy parameter of charged particles, v 2 , has 
been measured with the cumulant method using the ATLAS detector at the LHC. In a data sample 
corresponding to an integrated luminosity of approximately 1 /ib _1 , the parameter v 2 has been ob- 
tained using two- and four-particle cumulants over the pseudorapidity range \rj\ < 2.5. The results 
are presented as a function of transverse momentum and the event activity, defined in terms of the 
transverse energy summed over 3.1 < rj < 4.9 in the direction of the Pb beam. They show features 
characteristic of collective anisotropic flow, similar to that observed in Pb+Pb collisions. A comparison 
is made to results obtained using two-particle correlation methods, and to predictions from hydrody- 
namic models of p+Pb collisions. Despite the small transverse spatial extent of the p+Pb collision 
system, the large magnitude of v 2 and its similarity to hydrodynamic predictions provide additional 
evidence for the importance of final-state effects in p+Pb reactions. 
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to hydrodynamic predictions provide additional evidence for the importance of final-state effects in p+Pb 
reactions. 



1. Introduction n 

Recent observations of ridge-like structures in 

29 

the two-particle correlation functions measured in 
proton-lead (p+Pb) collisions at 5.02 TeV [TH2] 30 
have led to differing theoretical explanations. These 

32 

structures have been attributed either to mecha- 

33 

nisms that emphasize initial-state effects, such as 

34 

the saturation of parton distributions in the Pb- 
nucleus @HZ], or to final-state effects, such as jet- 
medium interactions [8], interactions induced by 
multiple partons [9Tll2|. and collective anisotropic 

flow mCEB]. 

The collective flow of particles produced in nu- 
clear collisions, which manifests itself as a sig- 
nificant anisotropy in the plane perpendicular to 
the beam direction, has been extensively studied 

in heavy-ion experiments at the LHC [T5H2~1] and « 

RHIC (for a review see Refs. [25J US])- In P+Pb « 

collisions the small size of the produced system « 

compared to the mean free path of the interacting 45 

constituents might have been expected to generate 46 

weaker collective flow, if any, compared to heavy- 47 

ion collisions. 4a 

However, two-particle correlation studies per- 49 

formed recently on data from p+Pb collisions at 50 

v /snn = 5.02 TeV revealed the presence of a 51 



"ridge" , a structure extended in the relative pseudo- 
rapidity, A77, while narrow in the relative azimuthal 
angle, A<f>, on both the near-side (Acf> ~ 0) [T] and 
away-side (Acf> ~ tt) [21 [3]. Furthermore, it was 
shown in Refs. (2J.EI that, after subtracting the com- 
ponent due to momentum conservation, the A<p dis- 
tribution in high-multiplicity interactions exhibits a 
predominantly cos(2A0) shape, resembling the el- 
liptic flow modulation of the A(f> distributions in 
Pb+Pb collisions. 

The final-state anisotropy is usually character- 
ized by the coefficients, v n , of a Fourier decompo- 
sition of the event-by-event azimuthal angle distri- 
bution of produced particles [551 HZ] : 

v n = (cosn(</> - *„)), (1) 

where 4> is the azimuthal angle of the particle, \I/ n 
is the event-plane angle for the n-th harmonic, and 
the outer brackets denote an average over charged 
particles in an event. In non-central heavy- ion col- 
lisions, the large and dominating v-i coefficient is 
associated mainly with the elliptic shape of the nu- 
clear overlap, and #2 defines the direction which 
nominally points in the direction of the classical 
impact parameter. In practice, initial-state fluctu- 
ations can blur the relationship between VE^ and 
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52 the impact parameter direction in nucleus-nucleus 90 

53 collisions. In contrast, '5 2 in proton- nucleus would 91 

54 be unrelated to the impact parameter and deter- 92 

55 mined by the initial-state fluctuations. In nucleus- 93 

56 nucleus collisions, the V2 coefficient in central col- 94 

57 lisions and the other v n coefficients in all collisions 95 
sa are related to various geometric configurations aris- 96 

59 ing from fluctuations of the nucleon positions in the 97 

60 overlap region [28]. 98 

61 In this Letter, a direct measurement of the 99 

62 second-order anisotropy parameter, v% % is presented ™ 

63 for p+Pb collisions at ^/snn = 5.02 TeV. The cu- 101 

64 mulant method [2"9Tl3"2"] is applied to derive using 102 

65 two- and four-particle cumulants. The cumulant 103 

66 method has been developed to characterize true 104 

67 multi-particle correlations related to the collective 105 
ea expansion of the system, while suppressing correla- 106 

69 tions from resonance decays, Bose-Einstein corre- 107 

70 lations and jet production. Emphasis is placed on 10s 

71 the estimate of V2, ^{4}, obtained from the four- 109 

72 particle cumulants which are expected to be free no 

73 from the effects of short-range two-particle correla- m 

74 tions, e.g. from resonance decays, unlike the two- 112 

75 particle cumulants, used to estimate f2{2}. 113 

76 The measurements of multi-particle cumulants 114 

77 presented in this Letter should provide further con- us 

78 straints on the origin of long-range correlations ob- 116 

79 served in p+Pb collisions. i" 

118 
119 

so 2. Event and track selections 120 

121 

si The p+Pb data sample was collected during a 122 

82 short run in September 2012, when the LHC deliv- 123 

83 ered p+Pb collisions at the nucleon-nucleon centre- 124 

84 of-mass energy ^/snn = 5.02 TeV with the centre- 125 

85 of-mass frame shifted by —0.47 in rapidity relative 126 

86 to the nominal ATLAS coordinate framcO 127 

87 The measurements were performed using the 128 
ss ATLAS detector [33 . The inner detector (ID) 129 
89 was used for measuring trajectories and momenta 130 

131 
132 

1 ATLAS uses a right-handed coordinate system with its 133 

origin at the nominal interaction point (IP) in the centre of 134 
the detector and the z-axis along the beam pipe. The ic-axis 
points from the IP to the centre of the LHC ring, and the 

j/-axis points upward. Cylindrical coordinates (r, <j>) are used 136 

in the transverse plane, being the azimuthal angle around 137 

the beam pipe. For the p+Pb collisions, the incident Pb 138 
beam travelled in the +z direction. The pseudorapidity is 
defined in laboratory coordinates in terms of the polar angle 

6 as r) = — lntan(0/2). Transverse momentum and energy 140 

are defined as pt = ps'mB and = EsinS, respectively. 141 



of charged particles for \r)\ < 2.5 with the sil- 
icon pixel detector and silicon microstrip detec- 
tors (SCT), and a transition radiation tracker, all 
placed in a 2 T axial magnetic field. For event trig- 
gering, two sets of Minimum Bias Trigger Scintil- 
lators (MBTS), located symmetrically in front of 
the endcap calorimeters, at z = ±3.6 m and cov- 
ering the pseudorapidity range 2.1 < \r)\ < 3.9, 
were used. The trigger used to select minimum- 
bias p+Pb collisions requires a signal in at least 
two MBTS counters. This trigger is fully effi- 
cient for events with more than four reconstructed 
tracks with p^ > 0.1 GeV. The forward calorime- 
ters (FCal), consisting of two symmetric systems 
with tungsten and copper absorbers and liquid ar- 
gon as the active material, cover 3.1 < \q\ < 4.9 and 
are used to characterize the overall event activity. 

The event selection follows the same require- 
ments as used in the recent two-particle correlation 
analysis [3]. Events are required to have a recon- 
structed vertex with its z position within ±150 mm 
of the nominal interaction point. Beam-gas and 
photonuclear interactions are suppressed by requir- 
ing at least one hit in a MBTS counter on each side 
of the interaction point and at most a 10 ns dif- 
ference between times measured on the two sides 
to eliminate through-going particles. To elimi- 
nate multiple p+Pb collisions (about 2% of collision 
events have more than one reconstructed vertex), 
the events with two reconstructed vertices that are 
separated in z by more than 15 mm are rejected. 
In addition, for the cumulant analysis presented 
here, it is required that the number of reconstructed 
tracks per event, passing the track selections as de- 
scribed below, is greater than three. With all the 
above selections, the analysed sample consists of 
about 1.9 x 10 6 events. 

Charged particle tracks are reconstructed in the 
ID using the standard algorithm optimized for p+p 
minimum-bias measurements |34j . Tracks are re- 
quired to have at least six hits in the SCT detector 
and at least one hit in the pixel detector. A hit in 
the first pixel layer is also required when the track 
crosses an active pixel module in that layer. Addi- 
tional requirements are imposed on the transverse 
(do) and longitudinal (zq sin 9) impact parameters 
measured with respect to the primary vertex. These 
are: \do\ and \z sin(9| must be smaller than 1.5 mm 
and must satisfy Ido/a^J < 3 and |zosin#/er z | < 3, 
where ad and a z are uncertainties on the trans- 
verse and longitudinal impact parameters, respec- 
tively, as obtained from the covariance matrix of 
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Fig. 1: Upper plot: the EiSp b distribution with the six ac- 
tivity intervals indicated. Lower plot: the distribution of 
7V5Jj c for each activity interval. The leftmost distribution 
corresponds to the interval with the lowest EE^ b , etc. 



142 the track fit. The analysis is restricted to charged 

143 particles with 0.3 < pt < 5.0 GeV and < 2.5. 

144 The tracking efficiency is evaluated using 
«5 HIJING-generated [35] p+Pb events that are fully 

146 simulated in the detector using GEANT4 [Ml 137] . 

147 and processed through the same reconstruction 
us software as the data. The efficiency for charged 

149 hadrons is found to depend only weakly on the 

150 event multiplicity and on p^ for transverse mo 

151 menta above 0.5 GeV. An efficiency of about 82% 

152 is observed at mid-rapidity, < 1, decreasing to 

153 about 68% at \r/\ > 2. For low-px tracks, between 

154 0.3 GeV and 0.5 GeV, the efficiency ranges from 

155 74% at r) = to about 50% for \rj\ > 2. The number 

156 of reconstructed charged particle tracks, not cor- 

157 rected for tracking efficiency, is denoted by N^ c . 

158 The analysis is performed in different intervals 

159 of £E^ b , the sum of transverse energy measured 

160 in the FCal with 3.1 < rj < 4.9 in the direction of 

161 the Pb beam with no correction for the difference in 
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Table 1: Characterization of activity intervals as selected 
by Ei£^ b . In the last column, the mean and RMS of the 
number of reconstructed charged particles with \r}\ < 2.5 and 
0.3 < Pt < 5 GeV, N^°, is given for each activity interval. 



162 response to electrons and hadrons. The distribution 

163 of Si?^ b for events passing all selection criteria is 

164 shown in Fig. [T] These events are divided into six 

165 Y,E^ h intervals to study the variation of Vi with 
we overall event activity, as indicated in Fig. [T] and 

167 shown in Table [l] Event "activity" is characterized 

168 by ££^ b : the most active events are those with 

169 the largest T,E^ h . The distribution of N^ c for each 

170 activity interval is shown in the lower plot of Fig. [T] 

171 3. Data analysis 

172 The cumulant method involves the calculation of 

173 2fc-particle azimuthal correlations, corr n {2k}, and 

174 cumulants, c„{2fc}, where k = 1,2 for the analysis 

175 presented in this paper. The two- and four-particle 

176 correlations are defined as corr n {2} = (e 1 ™^ 1- ^ 2 -') 

177 and corr n {A} = { e M4>x+^-4>3-4>4,)^ respectively, 

178 where the angle brackets denote the average in 

179 a single event over all pairs and all combina- 
lao tions of four particles. After averaging over 

181 events, the two-particle cumulant is obtained as 

182 c„{2} = (corr„{2}), and the four-particle cumulant 

183 c„{4} = (corr„{4}) - 2 • (corr n {2}) 2 . Thus the ef- 

184 feet of two-particle correlations is explicitly re- 

185 moved in the expression for c„{4}. Further details 

186 are given in Refs. [29 l l30 l 132] . 

187 Direct calculation of multi-particle correlations 
las requires multiple passes over the particles in an 

189 event, and requires extensive computing time in 

190 high-multiplicity events. To mitigate this, it has 

191 been proposed in Ref. |32j to express multi-particle 

192 correlations in terms of the moments of the flow 

193 vector Q n , defined as Q n — X)i el "^S where the in- 

194 dex n denotes the flow harmonic and the sum runs 



3 



over all particles in an event. This analysis is re- 
stricted to the second harmonic coefficient, n = 2. 
The method based on the flow-vector moments en- 
ables the calculation of multi-particle cumulants in 
a single pass over the full set of particles in each 
event. 

The cumulant method involves two main steps 
[2U1 [3D]. In the first step, the so-called "refer- 
ence" flow harmonic coefficients are calculated us- 
ing multi-particle cumulants for particles selected 
inclusively from a broad range in px an d fj as: 



,,re£ 



,,re£ 



{2} = v^{2}, 
{4} = ^=^{4}, 



«2{2}(pr,?7) 



d 2 {2} 



(4) 



(2) 
(3) 

206 where V2 Cf {2} (i>2° f {4}) denotes the reference es- 

207 timate of the second-order anisotropy parameter 

208 obtained using two-particle, C2{2} (four-particle, 

209 c 2 {4}) cumulants. 

210 The flow-vector method is easiest to apply when 
2n the detector acceptance is azimuthally uniform [32 . 

212 A correction for any azimuthal non-uniformity in 

213 the reconstruction of charged particle tracks is ob- 

214 tained from the data |25j . based on an rj-cj) map 

215 of all reconstructed tracks. For each small (5r) = 

216 0.1, 8(f) — 27r/64) bin (labelled i), a weight is cal- 
21? culated as Wi(r),(j)) — (N(8i]))/Ni(Sr],5cj)), where 

218 (N(6r))) is the event-averaged number of tracks 

219 in the Srj slice to which this bin belongs, while 

220 Ni(5ri,S4)) is the number of tracks in an event 

221 within this bin. Using this weight forces the az- 

222 imuthal angle distribution of reference particles to 

223 be uniform in (f>, but it does not change the 77 

224 distribution of reconstructed tracks. A weighted 

225 Q-vector is evaluated as Q n = J2i w i em ^ i [32] • 

226 From Eqs. ^ and Q it is clear that the cumu- 

227 lant method can be used to estimate i>2 only when 

228 C2{4} is negative and C2{2} positive. 

229 In the second step, the harmonic coefficients are 

230 determined as functions of px and 77, in bins in each 

231 variable (10 bins of equal width are used in 77 and 22 

232 bins of varied width in pr)- These differential flow 

233 harmonics are calculated for "particles of interest" 

234 which fall into these small bins. First, the differen- 

235 tial cumulants, (^{2} and c?2{4}, are obtained by 

236 correlating every particle of interest with one and 

237 three reference particles respectively. The differen- 

238 tial second harmonic, t>2{2fc}(px, where k — 1, 2, 

239 is then calculated with respect to the reference flow 

240 as derived in Refs. l29l 1301: 



250 
251 
252 
253 
254 
255 
256 
257 
258 
259 



«2{4}(pr,?7) 



-rf 2 {4} 



3/ </=^{4}' 



(5) 



The differential i>2 harmonic is then integrated 
over wider phase-space bins, with each small bin 
weighted by the appropriate charged particle mul- 
tiplicity. This is obtained from the reconstructed 
multiplicity by applying 77- and px-dependent effi- 
ciency factors, determined from Monte Carlo (MC) 
simulation as discussed in the previous section. Due 
to the small number of events in the data sample, 
the final results are integrated over the full accep- 
tance in 77. 
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Fig. 2: The two-particle (upper plot) and four-particle (lower 
plot) cumulants calculated using the reference flow particles 
as a function of T,E^, h for data (circles), the fully simulated 
HUING events (open squares) and the large generator-level 
HUING sample (filled squares). For clarity, the points for 
the fully simulated (generated) HUING events are slightly 
shifted to the left (right). 



Fig. [2] shows the two- and four-particle cumu- 
lants, averaged over events in each event-activity 
class defined in Table [l] as a function of T,E^ h . It 
is observed that four-particle cumulants are neg- 
ative only in a certain range of event activity. 
This restricts subsequent analysis to events with 
Y,E^ h > 25 GeV, for which the four-particle cumu- 
lant in data is found to be less than zero by at least 
two standard deviations (statistical errors only). It 



260 was also checked that for these events C2{4} is un- 312 

261 changed within errors for any high-multiplicity se- 313 

262 lection. For example, defining N20 as the value of 314 

263 V c r f such that 20% of events have N£ c < N 20 315 

264 (i.e. N 2 q is the 20th percentile), then selecting 316 

265 N^ c > N20 leaves C2{4} unchanged within errors. 317 

266 And for 'SE^ h > 25 GeV this holds for any per- 31a 

267 centile selection. 319 

268 Fig. [2] also shows the cumulants calculated for 320 

269 50 million HIJING-generated events, using the true 321 

270 particle information only, as well as for one million 322 

271 fully simulated and reconstructed HIJING events, 323 

272 using the same methods as used for the data. The 324 

273 S_E^ b obtained from the HIJING sample is rescaled 325 

274 to match that measured in the data. It should be 326 

275 noted that the HIJING Monte Carlo model does 327 

276 not contain any collective flow, and the only corre- 32a 

277 lations are those due to resonance decays, jet pro- 329 

278 duction and momentum conservation. The values 330 

279 of C2{2} for HIJING events are smaller than the 331 

280 values obtained from the data, and there is no sig- 332 

281 nificant difference between the HIJING results ob- 333 

282 tained at the generator ("truth") level and at the 334 

283 reconstruction level. For C2{4}, the HIJING events 

284 at J^E^. h ~ 20 GeV show a negative value compa- 33j 

285 rable to the values seen in the data, indicating that 

286 correlations from jets or momentum conservation 336 

287 contribute significantly to i>2{4} in events of low 337 

288 multiplicity. For Y,E^ h > 25 GeV the generator- 33a 

289 level HIJING sample's values for C2 {4} are also neg- 339 

290 ative, but the magnitude is much smaller than in 340 

291 the data or in HIJING events with smaller ££^ b . 341 

292 The size of the fully simulated HIJING event sam- 342 

293 pie is too small to draw a definite conclusion about 343 

294 the sign or magnitude of C2{4}. 344 

295 The systematic uncertainties on «2{2} and «2{4} 345 

296 as a function of px and ££^ b have been evaluated 346 

297 by varying several aspects of the analysis proce- 347 

298 dure. Azimuthal-angle sine terms in the Fourier 348 

299 expansion should be zero, but a non-zero contribu- 349 

300 tion can arise due to detector biases. It was found 350 

301 that the magnitude of the sine terms relative to 351 

302 the cosine terms is negligible (below 1%) for ^{2} 352 

303 measured as a function of px, as well as for the 353 

304 px-integrated «2{2} and t>2{4}. In the case of the 354 

305 measurement of the px-dependent ^2{4}, the sys- 355 

306 tematic uncertainty attributed to the residual sine 356 

307 terms varies between 6% and 14% in the different 357 

308 Y,E^ h intervals. Uncertainties related to the track- 35a 

309 ing are obtained from the differences between the 359 

310 main results and those using tracking requirements 360 

311 modified to be either more or less restrictive. They 36i 



are found to be negligible (below 0.2%) for t>2{2}. 
For the px-dependent t>2{4} they give a contribu- 
tion of less than 6% to the systematic uncertainty, 
and less than 1% for the px-integrated «2{4}. In ad- 
dition to varying the track quality requirements, an 
uncertainty on the px dependence of the efficiency 
corrections is also taken into account, and found to 
be below 1% for the «2{2} and t>2{4} measurements. 
The correction of the azimuthal-angle uniformity is 
checked by comparing the results to those obtained 
with all weights, u>j, set equal to one. This change 
leads to small relative differences, below 1%, for the 
«2{2} measured as a function of px, as well as for 
the px-integrated ^{2} and i>2{4}. Up to 4% dif- 
ferences are observed in the px-dependent t>2{4}. 
All individual contributions to the systematic un- 
certainty are added in quadrature and quoted as the 
total systematic uncertainty. The total systematic 
uncertainties are below 1% for the t>2{2} measure- 
ment. The V2{4} measurement precision is limited 
by large statistical errors, whereas the systematic 
uncertainties stay below 15% for V2{4}(px) and be- 
low 2% for the px-integrated t>2{4}. 

4. Results 

Fig. [3] shows the transverse momentum depen- 
dence of W2{2} and «2{4} in four different classes 
of the event activity, selected according to ££^ b . 
A significant second-order harmonic is observed. 
«2{4} is systematically smaller than «2{2}, con- 
sistent with the suppression of non-flow effects in 
«2{4}. This difference is most pronounced at high 
Pt and in collisions with low Y,E^° where jet-like 
correlations not diluted by the underlying event can 
contribute significantly. Thus, i>2{4} appears to 
provide a more reliable estimate of the second-order 
anisotropy parameter of collective flow. As a func- 
tion of transverse momentum the second-order har- 
monic, «2{4}, increases with px up to px ~ 2 GeV. 
Large statistical errors preclude a definite conclu- 
sion about the px dependence of «2{4} at higher 
transverse momenta. 

The shape and magnitude of the px-dependence 
of i>2{4} is found to be similar to that observed 
in p+Pb collisions using two-particle correlations 
[21 [3]. The second-order harmonic, 1)2, can be ex- 
tracted from two-particle azimuthal correlations us- 
ing charged particle pairs with a large pseudorapid- 
ity gap to suppress the short-range correlations on 
the near-side (Acf> ~ 0) (3] [22] . However, the two- 
particle correlation measured this way may still be 
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Fig. 3: The second-order harmonic calculated with the two-particle (circles) and four-particle (stars) cumulants as a function 
of transverse momentum in four different activity intervals. Bars denote statistical errors; systematic uncertainties are shown 
as shaded bands. The V2 derived from the Fourier decomposition of two-particle correlations [3] is shown by squares. 



362 affected by the dijet correlations on the away-side 

363 (A(/> ~ 7r), which can span a large range in Ar). 

364 In Ref. [3], the away-side non-flow correlation is 

365 estimated using the yield measured in the lowest 

366 S_E^ b collisions and is then subtracted from the 

367 higher Y,E^ h collisions. The result of that study, 

368 «2{2PC}, is shown in Fig. [3l for the four activ- 

369 ity intervals with largest YjE^! , and compared to 

370 «2{4}. Good agreement is observed between «2{4} 

371 and v 2 {2PC} for collisions with Y,E^ h > 55 GeV. 

372 For 'SE^ h < 55 GeV, the disagreement could be 

373 due either to the subtraction procedure used to ob- 

374 tain v 2 {2PC} or to non-flow effects in t>2{4}, or to 

375 a combination. 

376 The dependence on the collision activity of the 

377 second-order harmonic, integrated over 0.3 < pt < 

378 5 GeV, is shown in Fig. |4j The large magni- 

379 tude of w 2 {2} compared to t>2{4} suggests a sub- 

380 stantial contamination from non-how correlations. 

381 The value of t>2{4} is approximately 0.06, with lit— 

382 tie dependence on the overall event activity for 

383 T,E^ h >25 GeV. The extracted values of i> 2 {4} 

384 are also compared to the V2{2PC} values obtained 

385 from two-particle correlations. Good agreement is 

386 observed at large £E^ b , while at lower "EE^ h the 

387 V2{2PC} is smaller than t>2{4}, which may be due 

388 to different sensitivity of the two methods to non- 
389 flow contributions that become more important in 

390 low Y,E^ h collisions. Although t>2{4} is constructed 

391 to suppress local two-particle correlations, it may 

392 still include true multi-particle correlations from 

393 jets, which should account for a larger fraction of 

394 the correlated particle production in the events with 
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Fig. 4: The second-order harmonic, V2, integrated over px 
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atic uncertainties are shown as shaded bands. Also shown is 
V2{2PC} (squares) and predictions from the hydrodynamic 
model |18 | (triangles) for the same selection of charged par- 
ticles as in the data. 



395 the lowest If the HIJING results, shown in 

396 Fig. [2] were used to correct the measured cumulants 

397 for this non-flow contribution, the extracted «2{4} 

398 would be decreased by at most 10% for t>2{4} shown 

399 in Fig. [4] However, this correction is not applied to 

400 the final results. 

401 It is notable that the trend of the pt depen- 

402 dcncc of both w 2 {4} and w 2 {2PC} in p+Pb col- 

403 lisions resembles that observed for v-i measured 

404 with the event-plane method in Pb+Pb collisions 
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at snn = 2.76 TeV [HES], although with a mag- 
nitude between that observed in the most central 
and peripheral Pb+Pb collisions. While the trend 
is found to be nearly independent of the Pb+Pb 
collision geometry, the magnitude in Pb+Pb events 
depends on the initial shape of the colliding sys- 
tem, and has been modelled for p<r < 2 GeV using 
viscous hydrodynamics 

Harmonic flow coefficients in p+Pb collisions at 
^/snn = 5.02 TeV have also been predicted using 
viscous hydrodynamics, with similar initial condi- 
tions as the Pb+Pb calculations (TS]. The pre- 
dicted magnitude of the second-order harmonic]^] 
is compared to the measured i>2{4} and V2{2PC} 
in Fig. |4j It can be seen that the hydrodynamic 
predictions agree with our measurements over the 
S_E^ b range where the model predictions are avail- 
able. 

5. Conclusions 

ATLAS has measured the second harmonic coef- 
ficient in p+Pb collisions at ^/snn = 5.02 TeV us- 
ing two- and four-particle cumulants. A significant 
magnitude of vi is observed using both two- and 
four-particle cumulants, although «2{2} is consis- 
tently larger than «2{4}, indicating a sizeable con- 
tribution of non-flow correlations to «2{2}. The 
transverse momentum dependence of i>2{4} shows 
a behaviour similar to that measured in Pb+Pb 
collisions at */snn = 2.76 TeV. The magnitude of 
«2{4} increases with pt up to about 2-3 GeV. As 
a function of the collision activity, i<2{4} remains 
constant, at the level of about 0.06, for the colli- 
sions with YiE^ 3 >25 GeV, which corresponds to 
about 40% of the data. The measured «2{4} is 
found to be consistent with the second harmonic co- 
efficient extracted by the Fourier decomposition of 
the long-range two-particle correlation function for 
collisions with ££^ b > 55 GeV. Good agreement is 
also found with the predictions of a hydrodynamic 
calculation for p+Pb collisions. 

Extending previous results based only on two- 
particle correlations, the multi-particle cumulant 
results presented here provide additional evidence 



448 for the importance of final-state effects in the high- 

449 est multiplicity p+Pb reactions. Final-state effects 

450 may lead to collective flow similar to that observed 

451 in the hot, dense system created in high-energy 

452 heavy-ion collisions. 
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